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Microbial processes such as nitrification and anaerobic ammonium oxidation (anammox) 
are important for nitrogen cycling in marine sediments. Seasonal variations of archaeal 
and bacterial ammonia oxidizers (AOA and AOB) and anammox bacteria, as well as the 
environmental factors affecting these groups, are not well studied. We have examined the 
seasonal and depth distribution of the abundance and potential activity of these microbial 
groups in coastal marine sediments of the southern North Sea. This was achieved by 
quantifying specific intact polar lipids as well as the abundance and gene expression of their 
16S rRNA gene, the ammonia monooxygenase subunit A (amoA) gene of AOA and AOB, 
and the hydrazine synthase (hzsA) gene of anammox bacteria. AOA, AOB, and anammox 
bacteria were detected and transcriptionally active down to 12 cm sediment depth. In all 
seasons, the abundance of AOA was higher compared to the AOB abundance suggesting 
that AOA play a more dominant role in aerobic ammonia oxidation in these sediments. 
Anammox bacteria were abundant and active even in oxygenated and bioturbated parts 
of the sediment. The abundance of AOA and AOB was relatively stable with depth and 
over the seasonal cycle, while anammox bacteria abundance and transcriptional activity 
were highest in August. North Sea sediments thus seem to provide a common, stable, 
ecological niche for AOA, AOB, and anammox bacteria. 
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INTRODUCTION 

Nitrogen is an essential component for living organisms and 
thus plays a critical role in controlling primary production 
(Gruber and Galloway, 2008). Microbial aerobic and anaerobic 
nitrogen cycling processes such as nitrification, denitrification and 
anaerobic ammonium oxidation are tightly coupled in marine sed- 
iments where 50% of marine nitrogen removal occurs (Zhang 
etal., 2013, and references therein). The discovery of ammonia- 
oxidizing archaea (AOA; Konneke et al., 2005; Wuchter et al., 2006) 
belonging to the phylum Thaumarchaeota (Brochier-Armanet 
etal., 2008), and bacteria of the Planctomycetes phylum per- 
forming the anaerobic ammonia oxidation (anammox) reaction 
(Strous etal., 1999; Kuypers etal., 2003), has changed our per- 
spective on the marine nitrogen cycle and the involved microbial 
key- players over the last decade. Ammonia oxidation is the first 
and rate-limiting step of nitrification in which ammonia is oxi- 
dized to nitrite by AOA and ammonia oxidizing bacteria (AOB; 
see Prosser and Nicol, 2008 for a review). The anammox pro- 
cess involves the combination of ammonium with nitrite to 
form dinitrogen gas which is then released from the system 
(Kuenen, 2008). 

Several studies focusing on ammonia oxidation in marine sed- 
iments have suggested that AOA have a more prominent role in 
marine nitrification than AOB based on the higher abundance of 



AOA amoA gene copies (Beman and Francis, 2006; Mosier and 
Francis, 2008; Santoro etal, 2008; Abell etal, 2010). Proposed 
environmental controls shaping the distribution and activity of 
AOA in marine and estuarine surface sediments include the pres- 
ence of sulfide, salinity, dissolved oxygen (DO) concentration, the 
availability of phosphorous, temperature, and primary produc- 
tion (Mosier and Francis, 2008; Sahan and Muyzer, 2008; Erguder 
etal, 2009; Bernhard etal, 2010; Dang etal, 2010; Sakami, 2012; 
Bale etal, 2013). Beman etal. (2012) studied the seasonal dis- 
tribution of AOA and AOB amoA gene abundance in Catalina 
Harbor (CA, USA) down to 10 cm sediment depth and showed 
higher abundances of AOA and AOB amoA genes during summer 
compared to winter months as well as coinciding higher 15 NH4 + 
oxidation rates. For anammox bacteria, studies have shown that 
temperature, organic carbon content, nitrite concentration, water 
depth, sediment characteristics, and bioturbation play a role in 
their distribution, abundance, and activity in marine and estuar- 
ine sediments (Dalsgaard and Thamdrup, 2002; Engstrom etal, 
2005; Meyer etal, 2005; Dale etal, 2009; Jaeschke etal., 2009; 
Dang etal, 2010; Li etal, 2010; Neubacher etal, 2011, 2013; 
Laverock etal., 2013; Zhang etal, 2013). However, most of the 
studies have focused on the spatial distribution of anammox 
bacteria abundance and activity only in surface sediments (Dals- 
gaard and Thamdrup, 2002; Thamdrup and Dalsgaard, 2002; 
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Engstrom etal., 2005; Nicholls and Trimmer, 2009; Bale etal., 
2014). Limited studies have been focused on the abundance of 
anammox bacteria with sediment depth or the impact of sea- 
sonality on this group (Jaeschke etal., 2009; Zhang etal., 2013). 
Anammox bacteria abundance and activity has been previously 
detected using anammox specific ladderane biomarker lipids in 
combination with 15 N-labeling experiments in continental shelf 
and slope sediments of the Irish and Celtic Sea down to eight 
centimeters depth with spatial differences (Jaeschke etal., 2009). 
Recently, anammox bacteria abundance was shown to be higher 
during summer compared to the winter months in surface sed- 
iments of three hyper-nutrified estuarine tidal flats of Laizhou 
Bay (Bohai Sea, China) through the quantification of anammox 
bacteria 16S rRNA gene (Zhang etal, 2013). A study by Lave- 
rock etal. (2013) indicated that the temporal variation in the 
abundance of N- cycling functional genes is directly influenced 
by bioturbation activity varying between different bacterial and 
archaeal N-cycling genes. The relative abundance of AOB amoA 
genes is significantly affected by bioturbation whereas the relative 
abundance of AOA amoA seems to be controlled by other fac- 
tors (Laverock et al, 2013). Interestingly, bioturbation and mixing 
can extend the area of nitrate reduction in sediments leading to 
the production of nitrite which would fuel the anammox process 
(Meyer etal, 2005). 

The aim of our study was to determine differences in abundance 
and activity of ammonia oxidizers with sediment depth (down to 
12 cm depth) in marine sediments from the southern North Sea, 
a continental shelf sea. To address this issue, we quantified the 
abundance of 16S rRNA gene of ammonia oxidizer groups (AOA, 
AOB, and anammox bacteria), metabolic genes (amoA of AOA 
and AOB, and hzsA gene of anammox bacteria), and the intact 
polar membrane lipids specific for anammox bacteria (i.e., ladder- 
ane lipids, Sinninghe Damste etal, 2002b), and Thaumarchaeota 
(i.e., crenarchaeol; Sinninghe Damste etal., 2002a). The potential 
transcriptional activity of the target genes was analyzed to eluci- 
date the potential activity of the target organisms in the sediment. 
Our results suggest that both aerobic and anaerobic ammonia oxi- 
dizers coexist in the same niche and are actively involved in the 
nitrogen cycle in oxygenated and bioturbated sediments in the 
North Sea. 

MATERIALS AND METHODS 
STUDY SITE AND SAMPLING 

The Oyster ground, an almost circular depression with a max- 
imum depth of 50 m located in the southern North Sea 
(Weston etal., 2008), is a temporary deposition center for 
sediment (Raaphorst etal, 1998), which plays an important 
role in the carbon and nitrogen cycle in the region (Weston 
etal., 2008). The sedimentation of organic matter in this 
region leads to organic rich and muddy sediment, which sup- 
ports high levels of benthic fauna (Duineveld etal., 1991; 
van Raaphorst etal., 1992). 

Sediment cores were collected at a station (4°33.01 / E, 54° 13.00' 
N) in the Oyster ground during three cruises on board of the 
R/V Pelagia in February, May, and August 2011. Sediment was 
collected with 10 cm diameter multicores. Bottom water of the 
overlaying water was collected using a syringe and filtered (through 



a 0.45 |xm 25 mm Acrodisc HT Tuffryn Membrane syringe filter) 
for nutrient analysis. The cores were sliced into 1 cm slices using 
a hydraulic slicer or a manual slicer. Samples were collected for 
lipid and DNA/RNA analysis and kept at -80°C (DNA/RNA) and 
— 40° C (lipids) until processing. 

Bottom water temperature, water depth and salinity were mea- 
sured using a SBE911+ conductivity- temperature depth (CTD) 
system (Seabird Electronics, Inc., WA, USA). The oxygen concen- 
tration in the water column was measured using a DO sensor SBE 
43 (Seabird Electronics, Inc., WA, USA) integrated in the CTD 
system. 

PHYSICOCHEMICAL PARAMETERS 

Pore water was extracted from 2.5 mm sediment slices from 
0 to 2 cm depth, and from 1 cm sediment slices from 2 to 
12 cm depth by centrifugation (~4000 x g, 5 min, through a 
0.45 |xm 25 mm Acrodisc HT Tuffryn Membrane syringe fil- 
ter), and the obtained concentrations of the first 2 cm were 
averaged for the first and second centimeter of the sediment to 
estimate pore water nutrient concentrations in a 1 cm resolu- 
tion. Pore water samples were stored in pre-rinsed pony vials 
and NH4 + , NO3 - , NO2 - , and PO4 3- concentrations were 
analyzed as described by Bale etal. (2013). The total organic 
carbon (TOC) and total organic nitrogen (TON) content per 
gram of freeze dried sediment were determined after acidifica- 
tion of freeze dried sediment (0.5—1 g) with 2 N HC1. Residues 
were analyzed by using a Thermo Finnigan Delta plus isotope 
ratio monitoring mass spectrometer (irmMS) connected to a 
Flash 2000 elemental analyzer (Thermo Fisher Scientific, Milan, 
Italy). 

INTACT POLAR LIPID EXTRACTION AND ANALYSIS 

Intact polar lipids (IPLs) were extracted following Bale etal. 
(2013). Hexose-phosphohexose (HPH) crenarchaeol, a specific 
biomarker lipid for AOA (Schouten etal., 2008), was analyzed 
using high performance liquid chromatography mass spectrom- 
etry (HPLC/MS) as described by Pitcher etal. (2011a). Due 
to the lack of a standard the results are reported here as 
response unit per gram of sediment (r.u. g _1 ). The anam- 
mox bacteria specific intact ladderane phospholipid, C20- 
[3]-monoether ladderane, attached to a phosphatidylcholine 
(PC) headgroup (PC-monoether ladderane) was analyzed by 
HPLC/MS following Jaeschke etal. (2009) and quantified using 
an external standard consisting of isolated PC-monoether 
ladderane. 

DNA/RNA EXTRACTION 

DNA and RNA from sediment cores were extracted by using the 
DNA and RNA PowerSoil® Total Isolation Kit, respectively (Mo 
Bio Laboratories, Inc., Carlsbad, CA, USA). Nucleic acid con- 
centrations were quantified spectrophotometrically (Nanodrop, 
Thermo Scientific, Wilmington, DE, USA) and checked by agarose 
gel electrophoresis for integrity. Extracts were kept frozen at 
— 80° C. The RNA extracts were treated with RNase-free DNase 
(DNA- free™, Ambion Inc., Austin, TX, USA). RNA quality 
and concentration were estimated by the Experion RNA Std- 
Sens Analysis Kit (Bio-Rad Laboratories, Hercules, CA, USA). 
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DNA contamination was checked by PCR using RNA as a 
template. 

REVERSE TRANSCRIPTION (RT)-PCR 

Reverse transcription (RT) was performed with the Enhanced 
Avian First Strand synthesis kit (Sigma-Aldrich Co., St 
Louis, MO, USA) using random nonamers as described pre- 
viously (Holmes etal., 2004). Two negative controls lack- 
ing reverse transcriptase or RNA were included. PCR 
reactions were performed as described above to confirm 
the transcription to complementary DNA (cDNA) and 
the negative controls using the RT reaction as a tem- 
plate. 

QUANTITATIVE PCR (qPCR) ANALYSIS 

Quantitative PCR (qPCR) analyses were performed on a Biorad 
CFX96™ Real-Time System/ClOOO Thermal cycler equipped with 
CFX Manager™ Software. Gene copy numbers of the Thau- 
marchaeota group 1.1a were estimated by using the 16S rRNA 
gene primers MCGI-391F/MCGI-554R (Coolen etal, 2007). The 
AOA amoA gene was amplified using the primer combina- 
tion AmoA-ModF/AmoA-ModR (Yakimov etal, 2011). Abun- 
dance of AOB 16S rRNA gene was estimated using primers 
CT0189F/CT0654R as described by Kowalchuk etal. (1997). 
Gene copy numbers of the AOB amoA gene were estimated using 
the primer set AOB-araoAF/AOB-araoAR new (Rotthauwe etal., 
1997; Hornek etal., 2006). Abundance of the anammox bacteria 
16S rRNA gene was estimated using primers Brod541F/Amx820R 
as described by Li etal. (2010). Additionally, the anammox bac- 
teria hzsA gene was quantified using the primer combination 
/zzsA_1597F//zzsA_1857R as described by Harhangi etal. (2012; 
see Table SI in Supplementary Materials for details). All qPCR 
reactions were performed in triplicate with standard curves from 
10° to 10 7 molecules per microliter. Standard curves were gen- 
erated as described before (Pitcher etal., 2011b). Gene copies 
were determined in triplicates on diluted DNA extracts (1:10) 
and on cDNA extracts. The reaction mixture (25 |xL) contained 
1 U of PicoMaxx high fidelity DNA polymerase (Stratagene, 
Agilent Technologies, Santa Clara, CA, USA) 2.5 |xL of lOx 
PicoMaxx PCR buffer, 2.5 |xL 2.5 mM of each dNTP, 0.5 uX 
BSA (20 mg rnL -1 ), 0.02 pmol |xL _1 of primers, 10,000 times 
diluted SYBR Green® (Life Technologies, Carlsbad, CA, USA; 
optimized concentration), 0.5 |xL 50 mM of MgCi2, and ultra- 
pure sterile water. All reactions were performed in iCycler iQ™ 
96-well plates with optical tape (Bio-Rad, Hercules, CA, USA). 
Specificity of the reaction was tested with a gradient melting 
temperature assay. The cycling conditions for the qPCR reac- 
tion were the following: 95°C, 4 min; 40-45 x [95°C, 30 s; 
melting temperature (T m ), 40 s; 72°C, 30 s]; final exten- 
sion 80° C, 25 s. Specificity for qPCR reaction was tested on 
agarose gel electrophoresis and with a melting curve analysis (50- 
95°C; with a read every 0.5°C held for 1 s between each read) 
in order to identify unspecific PCR products such as primer 
dimers or fragments with unexpected fragment lengths. Melt- 
ing temperature, PCR efficiencies (£), and correlation coefficients 
for standard curves are listed in Table S2 in Supplementary 
Materials. 



PCR AMPLIFICATION AND CLONING 

Amplifications of the anammox bacteria 16S rRNA gene, AOA, 
and AOB amoA genes were performed with the primer pairs spec- 
ified above (Table SI in Supplementary Materials). PCR reaction 
mixture was the following (final concentration): Q-solution (PCR 
additive, Qiagen, Valencia, CA, USA) lx; PCR buffer lx; BSA 
(200 |xg mL" 1 ); dNTPs (20 |xM); primers (0.2 pmol |xL _1 ); 
MgCi2 (1.5 mM); 1.25 U Taq polymerase (Qiagen, Valencia, 
CA, USA). PCR conditions for these amplifications were the fol- 
lowing: 95°C, 5 min; 35 x [95°C, 1 min; T m , 1 min; 72°C, 
1 min]; final extension 72°C, 5 min. PCR products were gel 
purified (QIAquick gel purification kit, Qiagen, Valencia, CA, 
USA), cloned in the TOPO-TA cloning® kit (Life Technolo- 
gies, Carlsbad, CA, USA), and transformed in E. coll TOP 10 
cells following the manufacturer's recommendations. Recombi- 
nant plasmidic DNA was sequenced using M13R ^-CAG GAA 
ACA GCT ATG AC-3 r ) primer by Macrogen Inc. (Amsterdam, 
Netherlands). 

PHYL0GENETIC ANALYSIS 

Sequences were analyzed for the presence of chimeras 
using the Bellerophon tool at the GreenGenes website 
(http://greengenes.lbl.gov/). Sequences were aligned with Mega5 
software (Tamura etal., 2011) by using the alignment method 
ClustalW. Partial sequences of the archaeal amoA gene generated 
in this study were added to the amoA gene reference tree pro- 
vided in Pester et al. (2012) using the ARB Parsimony tool (Ludwig 
etal., 2004). Phylogenetic trees of AOB amoA and the anammox 
bacterial 16S rRNA genes were computed with the Neighbor- 
Joining method (Saitou and Nei, 1987) in the Mega5 software. 
The evolutionary distances were estimated using the Jukes-Cantor 
method (Jukes and Cantor, 1969) with a bootstrap test of 1000 
replicates. Sequences were deposited in NCBI with the follow- 
ing accession numbers: KJ807530-KJ807556 for partial sequences 
of the archaeal amoA gene, KJ807557-KJ807597 for partial bac- 
terial amoA gene sequences and KJ807598-KJ807609 for partial 
anammox bacterial 16S rRNA gene sequences. 

STATISTICAL ANALYSIS 

Spearman's rank order correlation coefficient (r s ) analysis was 
performed using the SigmaPlot (12.0) Exact Graphs and Data 
Analysis (Systat Inc., San Jose, CA, USA). Additional multivariate 
analysis did not show additional correlations of environmental 
data with gene and expression data therefore results were not 
included in this study. 

RESULTS 

PHYSICOCHEMICAL CONDITIONS 

Oxygen concentrations of the water column corresponding to the 
sampling time and location of sediment core sampling are shown 
in Figure SI in Supplementary Materials. In February, oxygen con- 
centrations decreased throughout the water column from 301.3 to 
299 |xM (300 |xM on average). In May, the oxygen concentration 
was slightly lower compared to February (282 |xM on average) 
and decreased between 17.8 and 22.8 m water depth from about 
287 to 278 |xM. Lowest oxygen concentrations were detected in 
August compared to February and May (230 |xM on average) 
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and decreased between 17.8 and 24.8 m water depth from 242 
to 217.25 |xM. 

Physicochemical parameters of the bottom water during the 
cruises of February (winter), May (spring), and August (sum- 
mer) are shown in Table S3 in Supplementary Materials. Bottom 
water temperatures ranged from 5 to 15.4°C. Salinity values 
were stable throughout the year ranging from 34.3 to 34.7 
practical salinity units (psu). Ammonia (NH4 + ) concentra- 
tions in the bottom water ranged between 1.6 and 3 |xM with 
a maximum concentration in February. Bottom water nitrite 
(NO2 - ) concentrations varied from 0.6 |xM in February to 
0.1 |xM in May and August. Highest nitrate (N03~) concentra- 
tions were detected in February (1 |xM) and lowest in August 
(0.5 |xM). 

Pore water concentrations of NH4 + ranged between 3 and 
53 |xM in all three seasons (Figure 1). The depth trend of 
ammonia concentrations in the pore water was similar in all 
seasons with lower concentrations in the upper 4 cm below 
sea floor (bsf; 3-34 |xM) compared to increasing concentra- 
tions of the underlying layers (20-53 |xM) down to 12 cm 
bsf. Highest NH4 + concentrations were detected in August 
(40-53 |xM) compared to February and May. Nitrite (NC>2 _ ) 
concentrations ranged between 0.1 and 1.2 |xM in all seasons 
(Figure 1). Depth profiles of pore water nitrite concentrations 
showed different seasonal trends. In February, nitrite concen- 
trations were relatively low (0.3 |xM) in the first 3 cm bsf and 
increased with depth (0.6 |xM) whereas in May and August, 
nitrite concentrations were higher between 1 and 2 cm bsf 
(0.8 |xM) compared to the underlying sediment layers (0.4 and 
0.2 |xM, respectively). Pore water nitrate (NO3 - ) concentra- 
tions strongly varied with sediment depth and season (Figure 1). 
In February and May, nitrate concentrations were highly vari- 
able with depth and reached maximum values between 3 and 



4 cm (38 |xM) and 4 and 5 cm bsf (24 |xM), while the concen- 
tration strongly decreased in the underlying layers. In August, 
lower nitrate concentrations were detected compared to the 
other seasons with slightly elevated values in the first 2 cm 
bsf (7.3 |xM) and lower concentrations in the underlying lay- 
ers (1.3 |xM). Phosphate concentrations ranged between 1.4 and 
5.2 |xM in all seasons (Figure 1). Pore water phosphate con- 
centrations were stable in the first 4 cm bsf (1.7 |xM, averaged) 
in all three seasons. In February and May, phosphate concen- 
trations increased slightly (3.4 |xM, averaged) in the underlying 
sediment layers whereas the depth profile remained stable in 
August. TOC and TON percentages of the sediment were rel- 
atively stable in all seasons (between 0.15 and 0.38% and 0.02 
and 0.08%, respectively), with highest values in May between 
1 and 2 cm and 4 and 5 cm bsf (Table S4 in Supplementary 
Materials). 

ABUNDANCE, DISTRIBUTION, ACTIVITY, AND DIVERSITY OF A0A AND 
A0B 

The abundance of AOA in the sediments was determined by quan- 
tification of the 16S rRNA and amoA gene copies, as well as by 
quantification of HPH-crenarchaeol (Figures 2A-C). The sea- 
sonal variations in depth profiles of AOA 16S rRNA and amoA 
gene abundance were generally similar (Figures 2A,B). Both AOA 
gene abundances were higher between 0 and 5 cm (6 x 10 6 and 
1.4 x 10 7 gene copies g -1 , respectively) compared to 5 and 
12 cm bsf (5 x 10 6 and 8.2 x 10 6 gene copies g -1 , respec- 
tively). Gene copy numbers of the AOA amoA gene were higher 
in August in comparison to the values observed in February and 
May (1.4-fold and 2.7-fold higher, respectively) for the 0-5 cm 
interval. HPH-crenarchaeol concentration was variable over the 
first 5 cm bsf (between 3 x 10 6 and 17 x 10 6 r.u. g -1 ), but 
stable values (around 5.9 x 10 6 r.u. g _1 ) were detected in the 
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FIGURE 1 | Sediment pore water nutrients: (A) NH 4 +; (B) N0 2 ~; (C) N0 3 _ and (D) P0 4 3+ (jiM) with sediment depth (cm) below sea floor (bsf). 
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FIGURE 2 | Abundance of (A) ammonia oxidizing Archaea (AOA) 16S 
rRNA gene (copy number g~ 1 ); (B) AOA amok gene (copy number 
g~ 1 ); (C) Hexose phosphohexose (HPH)-crenarchaeol (response 



units g" 1 ), and (D) AOB 16S rRNA gene (copy number g" 1 ); (E) AOB 
amok gene (copy number g" 1 ) with sediment depth (cm) below sea 
floor (bsf). 



layers underneath. No clear seasonal differences were detected 
(Figure 2C). 

The abundance of AOB was determined by quantifying the 
AOB 16S rRNA and amoA gene (Figures 2D,E). Values were rel- 
atively stable with depth with seasonal variations between 0 and 
4 cm bsf. Abundance of the AOB 16S rRNA gene showed val- 
ues similar to the AOA 16S rRNA gene abundance in February 
and August (on average 5.1 x 10 6 and 5.8 x 10 6 gene copies 
g _1 , respectively), while in May AOB 16S rRNA gene copy num- 
bers were lower on average compared to the AOA 16S rRNA 
gene copy number (4.0 x 10 6 gene copies g _1 , 1.4-fold lower). 
The AOB amoA gene abundance followed the same seasonal 
trends as that of the AOB 16S rRNA gene, AOA 16S rRNA and 
amoA gene abundances, with higher values in August, but abso- 
lute values were one order of magnitude lower than AOA amoA 
gene copy numbers (2.3 x 10 5 -5.2 x 10 6 gene copies g _1 ). 
AOB amoA gene abundance was highest in the 0-5 cm depth 
interval except in May when values between 2-5 cm and in 
the 5-12 cm depth interval were lower than in the surface 
sediment. 

The RNA:DNA ratios of the AOA and AOB 16S rRNA and amoA 
genes were calculated as an indicator of the potential transcrip- 
tional activity of the targeted microbial group (Figures 3A-D). 
In all three seasons, AOA 16S rRNA gene RNA:DNA ratio 
showed relatively stable values down core ranging between 0.9 
and 5.0 and 1.2 and 3.4 in May and August but reached slightly 
higher values (2.2-8.2) in February (Figure 3A). AOB 16S rRNA 
gene RNA:DNA ratio was between 0.1 and 2.1 with a min- 
imum of 0.1 between 4 and 5 cm bsf in February and 0.1 
and 3.7 with a minimum between 2 and 3 cm bsf in May. 
The AOB 16S rRNA gene RNA:DNA ratio was slightly higher 



(1.2-4.1) in August (Figure 3C). Values were relatively sta- 
ble with depth in August whereas in February and May values 
were more variable. Compared to the RNA:DNA ratio of the 
AOA 16S rRNA gene, the AOB 16S rRNA gene RNA:DNA ratio 
was lower (fivefold in February, 1.5-fold in May and 1.3-fold 
lower in August) in all seasons (Figures 3A,C). AOA amoA gene 
transcripts (Figure 3B) were only detected in February in the 
first 3 cm and in August throughout the core, and were under 
the detection limit in May. AOB amoA gene RNA:DNA ratio 
(Figure 3D) varied between 4 x 10~ 5 and 2 x 10 -2 in all 
seasons with a slight decrease with depth and highest values in 
August. 

The diversity of AOA and AOB was evaluated in two selected 
depth intervals (0-1 and 9-10 cm bsf) of the sediment core recov- 
ered in August by targeting the AOA amoA gene and the AOB amoA 
gene, respectively (Figures S2 and S3 in Supplementary Materi- 
als). Amplified AOA amoA gene sequences were closely related to 
AOA amoA gene sequences previously isolated from marine and 
estuarine environments and did not cluster according to depth 
(Figure S2 in Supplementary Materials). Sequences were mainly 
affiliated to the Nitrosopumilus subclusters 12 and 16, also known 
as the stable marine cluster, or to the Nitrosopumilus subcluster 
4.1, known as the estuarine cluster, according to the phylogenetic 
classification proposed by Pester etal. (2012). The AOB amoA 
gene sequences were affiliated to AOB amoA gene sequences pre- 
viously isolated from estuarine sediments. Sequences were closely 
related to the AOB amoA gene sequences of Nitrospira sp. (acces- 
sion number X90821.1) and Nitrosolobus multiformis (accession 
number AF 042171.1). No clustering of AOB amoA sequences 
according to depth was observed (Figure S3 in Supplementary 
Materials). 
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FIGURE 3 | Transcriptional activity given as RNA:DNA ratio of (A) AOA 16S rRNA gene; (B) AOA amoA gene; (C) AOB 16S rRNA gene, (D) AOB amoA; 
(E) Anammox 16S rRNA gene, and (F) Anammox bacteria hzsA gene with sediment depth (cm) below sea floor (bsf). 



ABUNDANCE, DISTRIBUTION, ACTIVITY, AND DIVERSITY OF 
ANAMMOX BACTERIA 

The abundance and distribution of anammox bacteria were esti- 
mated by quantification of anammox bacteria 16S rRNA and 
hzsA gene copy number, as well as by quantification of the 
PC-monoether ladderane (Figures 4A-C). Depth profiles of 
anammox bacteria 16S rRNA and hzsA gene abundances fol- 
lowed similar trends. Gene copy numbers of the hzsA gene 
were one order of magnitude lower (between 3.5 x 10 5 and 
2 x 10 6 copies g _1 ) than anammox bacteria 16S rRNA gene 
copy numbers in all three seasons (Figures 4A,B). Values 
were higher in August compared to those in February and 
May, especially for the anammox bacteria 16S rRNA gene 
(3. 7 -fold higher). The RNA:DNA ratio for the anammox bac- 
teria 16S rRNA gene varied between 1.6 and 34.6 (Figure 3E) 
with higher values in August compared to February and May 
(on average fourfold higher). The hzsA gene RNA:DNA ratio 
(Figure 3F) depth profiles showed relatively stable values (between 
0.008 and 0.125) throughout the core without significant sea- 
sonal differences. The PC-monoether ladderane concentration 
(Figure 4C) was variable over the first 5 cm bsf (between 0.25 
and 1.25 ng g _1 ), while in the underlying sediment layers 
values decreased slightly to more stable values (approximately 
0.36 ng g -1 ) down to 12 cm bsf. There was no clear seasonal 
difference observed. 

The diversity of anammox bacteria was evaluated in two 
selected depths (0-1 and 9-10 cm bsf) of the sediment core in 
August by targeting the 16S rRNA gene of anammox bacteria. The 
sequences were closely related to "Candidatus Scalindula marina" 
(accession number EF602038) and "Candidatus Scalindula bro- 
dae" (accession number AY254883) and to sequences previously 
detected in waters of marine oxygen minimum zones (Woebken 



etal., 2008). Anammox bacterial 16S rRNA gene sequences from 
different depths were closely related to each other (Figure S4 in 
Supplementary Materials). 

DISCUSSION 

AEROBIC AMMONIA OXIDIZERS 

In this study, we defined two niches in the Oyster ground sedi- 
ment based on the differences in physicochemical parameters and 
in the abundance of the target organisms, i.e., the upper 4 cm of 
the sediment, and the underlying sediment. Ammonia concentra- 
tion in the pore water of the upper part of the sediment (0-4 cm) 
was consistently lower than in deeper sediment layers (4-12 cm; 
Figure 1), suggesting ammonia diffusion from below formed by 
mineralization of organic matter by ammonification or dissimi- 
latory nitrate reduction to ammonium (DNRA). In addition, the 
abundance of AOA and AOB (as well as anammox bacteria, see 
below) was higher and more affected by seasonality in the upper 
part than deeper in the sediment, which apparently is a more sta- 
ble compartment. The higher variability observed in the upper 
part of the sediment is in agreement with the seasonal mixing of 
the bottom water layers and sediment deposition dynamics in this 
area (Greenwood etal., 2010). 

Oxygen availability is expected to play an important role in 
the dynamics of ammonia oxidizing microorganisms in marine 
sediments. In the summer months, water stratification leads to 
a decrease in the DO concentration in the bottom water in the 
Oyster ground area (Greenwood et al., 2010; see Figure SI in Sup- 
plementary Materials). Indeed, bottom water hypoxia has been 
shown to cause a decrease in oxygen penetration depth and oxygen 
consumption in Oyster ground sediments by induced short-term 
hypoxia in intact sediment cores (Neubacher etal., 2011). Lohse 
etal. (1996) reported an oxygen penetration depth (OPD) of 
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FIGURE 4 | Abundance of (A) Anammox bacteria 16S rRNA gene (copy number g 1 ); (B) Anammox bacteria hzsA gene (copy number g 1 ), and (C) 
PC-monoether ladderane (ng g 1 ) with sediment depth (cm) below sea floor (bsf). 



1.7 cm in February and 0.5 cm in August in sediment at a station 
close to the Oyster ground using oxygen microelectrodes, suggest- 
ing permanently anoxic conditions from a depth of approximately 
>2 cm throughout the year. However, the OPD into the sediment 
must be interpreted with caution because the Oyster ground sedi- 
ment is known to be an area of intense bioturbation by macro- and 
meiobenthos even in periods of low oxygen concentration in the 
bottom water (De Wilde etal, 1984; Duineveld etal, 1991), thus 
the import of oxygen by macro faunal burrows into the deeper sed- 
iment could be relevant (Kristensen, 2000; Meysman etal., 2006). 
A study by Laverock etal. (2010) has also shown that bioturba- 
tion can induce changes in the microbial community composition, 
resulting in distinct communities within the burrow. However, in 
our study we did not observe changes in the diversity of amoA 
gene sequences of AOA or AOB, which might indicate that the 
changing physicochemical conditions in this system do not select 
for specific ammonia oxidizing microorganisms. The depth and 
seasonal variability for AOB was more pronounced than for AOA 
in the first 4 cm of the sediment, which is consistent with the fact 
that the relative abundance of AOB is more significantly affected 
by bioturbation than AOA according to previous studies (Laverock 
etal, 2013). 

AOA and AOB 16S rRNA and amoA genes and AOA-IPLs and 
were detected up to 12 cm depth in Oyster ground sediments. 
The detection of aerobic ammonia oxidizing microorganisms in 
marine sediment at depths where the oxygen concentration is 
expected to be undetectable seems counterintuitive with their aer- 
obic metabolism. However, the intense bioturbation reported in 
the Oyster ground sediments, as mentioned above, might provide 
enough oxygen to support the activity of these microbial groups. 



In fact, Beman etal. (2012) detected AOB and AOA amoA genes 
up to 10 cm sediment depth as well as readily detectable 15 NH4 + 
oxidation rates up to 9 cm bsf with maxima overlapping with a 
peak in pore water NO3 - + NO2 - concentrations between 3 and 
7 cm bsf in bioturbated sediments of Catalina Island (CA, USA). 
Thus, our results, and those of Beman etal. (2012), suggest that 
aerobic ammonia oxidizers are indeed living and actively involved 
in the marine sedimentary nitrogen cycle deeper (i.e., up to at least 
ca. 10 cm) in coastal marine sediments. 

The detection of a potential transcriptional activity of the 
AOA 16S rRNA gene up to 12 cm depth supports that AOA are 
indeed active deeper in the marine sediment. The potential AOA 
16S rRNA gene transcriptional activity was higher during win- 
ter, whereas the abundance was slightly higher during summer 
months. Previous studies showed higher AOA abundance during 
winter months in the in the North Sea water column (Wuchter 
etal., 2006; Herfort etal, 2007; Pitcher etal, 2011c) attributed 
to the higher availability of ammonia and the lack of competi- 
tion for ammonia with phytoplankton. However, none of these 
studies targeted the transcriptional activity of AOA. Recently, 
Bale etal. (2013) observed that the abundance of the benthic 
Thaumarchaeota (given by HPH crenarchaeol) in the surface sed- 
iment (0-1 cm) at the Oyster grounds was on average highest 
in the summer (August) and lowest in the winter (November), 
which they attributed to the deposition of algal-bloom-derived 
organic matter taking place in late spring and summer result- 
ing in the formation of ammonia. However, it should also be 
noted that Bale etal. (2013) did find that transcriptional activ- 
ity of the AOA 16S rRNA gene was higher during winter in 
the surface sediment. The disparity in seasonality between the 
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results of these two studies highlights the difference between 
the seasonal ammonia dynamics at the surface sediment/water 
interface (Bale etal., 2013) and with depth in the sediment (this 
study). 

To complement the determination of potential transcriptional 
activity of AOA based on 16S rRNA gene, we also determined the 
AO A amoA gene RNA:DNA ratio, which was possible for samples 
from August throughout the core but only in the upper 3 cm depth 
in February because the quantity of transcripts in the other seasons 
and depths was under the detection limit of the qPCR assay. 

The relationship of amoA gene expression and in situ ammo- 
nia oxidation is still unclear (Lam etal., 2007; Nicol etal., 2008; 
Abell et al., 201 1; Vissers et al, 2013). However, recent studies have 
observed a good correlation between nitrification rates and abun- 
dance of AOA amoA gene and transcripts in coastal waters (Smith 
et al, 2014; Urakawa et al, 2014). Therefore, the detection of AOA 
amoA transcripts in our setting suggests that AOA were active in 
the summer throughout the sediment in comparison with the win- 
ter in which AOA would be only active in the upper part of the 
sediment (upper 2 cm). 

AOB 16S rRNA gene abundance was in the same order of mag- 
nitude than AOA 16S rRNA gene abundance; however, AOB amoA 
gene abundance was one order of magnitude lower than AOA 
amoA gene abundance, which can be explained by biases intro- 
duced by the primers or by differences in copy number of 16S 
rRNA and amoA genes per cell. In fact, it is assumed that 3.3 ±1.6 
copies of 16S rRNA can exist per cell in AOB (Vetrovsky and 
Baldrian, 2013), whereas marine Crenarchaeota seem to contain 
only one copy per genome (Klappenbach etal, 2001; Laverock 
etal., 2013). For the amoA gene it is assumed that AOB genomes 
contain two to three copies (Arp et al., 2007) while AOA genomes 
contain only one copy. In addition, the transcriptional activity of 
the AOB 16S rRNA gene was lower compared to the one of AOA 
16S rRNA gene in all seasons. 

On the other hand, the abundance of AOB amoA gene tran- 
scripts was higher than for AOA and was detected up to 12 cm 
depth in all seasons which would suggest a higher AOB activity vs. 
AOA. However, a recent study by Urakawa etal. (2014) in water 
of the Puget Sound Estuary did not observe a correlation between 
nitrification rates and AOB amoA transcripts. In addition, previ- 
ous studies have detected the presence of AOB amoA transcripts in 
non-active starved cells (Bollmann etal., 2005; Geets etal., 2006; 
Sayavedra-Soto and Arp, 2011) indicating that the AOB amoA 
transcriptional activity is not a good indicator of AOB nitrifica- 
tion activity due to a longer half-life of AOB amoA transcripts than 
that expected for AOA amoA transcripts. 

Taking together the AOA and AOB amoA gene quantification 
results, AOA have a more important role in the nitrification process 
in Oyster ground sediments than AOB. In addition, all available 
evidence points to the presence of some oxygen deeper in the 
sediments of the Oyster ground in all seasons probably due to the 
intense bioturbation in the area, which favors the activity of both 
AOB and AOA. In the case of AOA, and due to their higher affinity 
for oxygen ( Martens -Habbena etal., 2009), the niche occupancy 
of AOA could be potentially larger than for AOB. 

Previous studies observed a positive correlation between AOA 
abundance and low phosphate concentrations and suggested that 



phosphate plays a role in determining the niches of Thaumar- 
chaeota (Erguder etal, 2009; Dang etal, 2013). In our study, 
AOA amoA gene transcript copy number and RNA:DNA ratio 
was negatively correlated with phosphate (r s = —0.745 and 
—0.739, respectively; both P < 0.005) suggesting a lower expres- 
sion of the AOA amoA gene with an increase in pore water 
phosphate concentrations. These results support the assumption 
that some ecotypes of AOA dominate over AOB in environ- 
ments with low phosphate availability. Previously, homologs of 
the pst gene, encoding for the Pst phosphate ABC transporter, 
which is usually activated in phosphate deficient environments, 
was found in marine metagenomes of thaumarchaeotal origin 
(Rusch et al., 2007; Dang et al., 2013) emphasizing the importance 
of phosphate availability in marine environments. No correla- 
tions were found between AOB abundance or transcriptional 
activity with pore water nitrogen species, but a negative corre- 
lation with the sediment depth (r s = -0.646; P < 0.005) was 
observed indicating that sediment depth is a factor determin- 
ing the AOB abundance and activity in Oyster ground sediments 
which might be explained by the decreasing oxygen availability 
with depth. 

ANAEROBIC AMMONIA OXIDIZERS 

Anammox bacteria 16S rRNA and hzsA gene abundance was 
detected throughout the analyzed sediment with more pro- 
nounced depth and seasonal differences than those observed for 
AOA and AOB. This was not reflected in the abundance of PC- 
monoether ladderane lipid. It is possible that the PC-monoether 
ladderane is partly derived from fossil material and thus reflects 
a long-term presence of anammox bacteria averaged throughout 
the entire year rather than an actively living anammox popula- 
tion. This possibility has been already suggested by Brandsma et al. 
(201 1) and Bale et al. (2014) who found dissimilarities between the 
PC-monoether ladderane concentration and gene copy numbers 
in surface sediments of the Gullmar Fjord and the Southern North 
Sea, respectively. 

We also observed an up to 10-fold difference in the quantifi- 
cation of anammox bacteria 16S rRNA and hzsA gene abundance 
(Figures 4A,B). This discordance between the anammox bacteria 
gene markers has been previously observed in environmental sam- 
ples (Harhangi etal, 2012), and it could be due to primer biases 
as previously suggested (Bale etal., 2014). Both the anammox 16S 
rRNA and the hzsA gene revealed a higher abundance in August 
compared to February and May. This may indicate that anam- 
mox bacteria are more active and abundant at higher temperatures 
(15°C), which is the anammox bacteria temperature optimum in 
temperate shelf sediments (Dalsgaard etal, 2005). Additionally, 
the oxygen penetration depth is expected to decrease during water 
stratification between May and October, which would also favor 
anammox bacterial growth and activity (Dalsgaard and Tham- 
drup, 2002). Previous studies by Neubacher et al. (2013) observed 
an increase in anammox activity up to 82% compared to con- 
trol treatments in mesocosms under sustained hypoxic conditions, 
which corresponds well with our findings of high anammox activ- 
ity during summer stratification associated with lower oxygen 
concentrations in the bottom water (Figure SI in Supplementary 
Materials). 
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The presence of anammox markers in the upper layers of the 
sediment (0-4 cm) is remarkable as the metabolism of anammox 
bacteria is expected to be incompatible with the presence of oxy- 
gen. However, some studies have reported that anammox bacteria 
can cope with low oxygen concentrations (Kalvelage etal, 2011; 
Yan etal, 2012), and that anammox bacteria can be present and 
active at high oxygen levels by being dormant or using anaerobic 
niches (Kuypers et al., 2005; Woebken et al., 2007). It is also possi- 
ble that a high activity of AOA and AOB and heterotrophic bacteria 
inhabiting these marine sediments would rapidly consume the 
oxygen available and allow the activity of anammox bacteria even 
in presence of oxygen. Some studies have also concluded that the 
potential for anammox bacteria activity is indeed not affected 
by bioturbation or physical mixing (Rooks etal, 2012). On the 
contrary, bioturbation and mixing can extend the area of nitrate 
reduction and thus the availability of nitrite, which fuels the anam- 
mox process (Meyer etal, 2005). Our results thus suggest that 
anammox bacteria can tolerate the presence of oxygen and co-exist 
with aerobic ammonia oxidizers and potentially have an important 
role in the nitrogen cycle in marine sediments. In fact, the poten- 
tial transcriptional activity of anammox bacteria 16S rRNA and 
hzsA gene was observed throughout the sediment in all seasons. 
The transcriptional activity of anammox bacteria 16S rRNA gene 
was higher in the summer as observed for the anammox bacteria 
abundance. A recent study by Bale etal. (2014) observed a good 
correlation between the rate of anammox and anammox bacte- 
ria 16S rRNA gene transcript abundance in North Sea sediments, 
which suggests that the transcriptional activity of anammox bac- 
teria 16S rRNA gene is a good proxy of anammox bacteria activity, 
and confirm a higher anammox activity in the summer for the sedi- 
ments analyzed in our study. Nevertheless, quantification results of 
anammox bacteria 16S rRNA gene transcripts have to be interpret 
with caution because there is evidence that the ribosome con- 
tent does not decrease during the period of starvation (Schmid 
etal., 2005; Li and Gu, 2011). However, the increase of anammox 
bacteria 16S rRNA gene transcriptional activity during summer 
was not observed for the hzsA gene, which remained stable across 
sediment depth and also for the different seasons, and gene abun- 
dances were generally an order of magnitude lower. The low and 
constitutive transcriptional activity of the hzsA gene indicates that 
the hzsA gene transcript abundance does not reflect variations 
of anammox bacteria activity in environmental sediment sam- 
ples, thus the hzsA gene transcript abundance does not seem 
to be an adequate biomarker for the estimation of the activity 
of anammox bacteria. Further experiments, especially with pure 
cultures and under controlled physiological conditions, should 
be performed to clarify the regulation of the expression of the 
hzsA gene. 

No evidence of a niche separation of aerobic and anaerobic 
ammonia oxidizers was observed in the oxygen transition zone 
of the marine sediments. Surprisingly, aerobic and anaerobic 
ammonia oxidizers are present and active throughout the year 
in oxygenated sediment layers as well as in anoxic layers. These 
results seem to be counterintuitive with regards to the individual 
oxygen requirements of the targeted microbial groups. We hypoth- 
esize that the presence and activity of aerobic ammonia oxidizers 
(AOA, AOB) is supported by the oxygen supply in deeper anoxic 



layers due to bioturbation. Likewise, the presence and activity of 
anammox bacteria in these sediments is then possible by the rapid 
consumption of oxygen by aerobic ammonia oxidizers and/or 
other groups, as well as the availability of nitrite provided either 
by ammonia oxidizers or by an active nitrate reduction present in 
bioturbated sediments (Meyer etal., 2005). 

CONCLUSION 

Our study has unraveled the coexistence and metabolic activ- 
ity of AOA, AOB, and anammox bacteria in bioturbated marine 
sediments of the North Sea, leading to the conclusion that 
the metabolism of these microbial groups is spatially cou- 
pled based on the rapid consumption of oxygen that allows 
the coexistence of aerobic and anaerobic ammonia oxidizers. 
AOA outnumbered AOB throughout the year which may be 
caused by the higher oxygen affinity of AOA compared to 
AOB. Anammox bacterial abundance and activity were higher 
during summer, indicating that their growth and activity are 
favored by higher temperatures and lower oxygen available in 
the sediments due to summer stratifying conditions in the 
water column. During the summer, anammox bacteria are 
probably not in competition with AOA and AOB for ammo- 
nia as the concentrations were relatively high in the sediment 
pore water most likely as a result of ammonification processes 
and the activity of denitrifiers and DNRA. Further studies 
are required to get a complete picture of the nature of these 
interactions in the oxygen transition zone of coastal marine 
sediments. 
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